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the  radiator  Is  quite  large.  Since  all  the  heat  must  be  radiated 
to  space,  the  radiator  also  gets  very  big  and  heavy. 

The  final  method  presently  available  for  transporting  waste 
heat  is  the  heat  pipe.  A  heat  pipe  is  nothing  more  than  a  tube, 
capped  at  both  ends,  that  contains  a  small  amount  of  liquid  and  a 
wick.  When  heat  is  added  to  one  end  (the  evaporator)  the  liquid 
boils.  Because  of  the  increase  in  vapor  pressure,  the  vapor  is 
transported  to  the  other  end  (the  condenser)  where  it  condenses  as 
heat  is  taken  out.  The  fluid  is  then  wicked  back  to  the  evaporator 
to  complete  the  cycle.  The  heat  pipe  is  an  extremely  efficient  and 
reliable  method  of  transporting  heat.  The  problem  with  heat  pipes, 
however,  is  that  they  can  only  transport  a  limited  amount  of  heat 
for  a  limited  distance^1".  It  is  possible  to  put  together  many 
heat-  pipe  bundles  in  series  to  overcome  this  problem,  but,  again 
weight  becomes  a  problem^. 

Due  to  the  shortcomings  of  the  previous  methods,  it  has  been 
proposed  to  use  two-phase-fluid  loops^.  Like  a  one-phase-fluid 
loop,  these  systems  would  involve  pumping  a  fluid  around  a  loop. 
Unlike  one-phase  systems,  two-phase  systems  would  use  the  latent 
heat  of  vaporization  rather  than  the  sensible  heat  of  the 
circulating  fluid.  Since  the  latent  heat  is  much  greater  than  the 
sensible  heat,  *Ms  system  can  use  a  much  smaller  flow  rate.  In 
fact,  preliminary  calculations  have  shown  that  the  pump  power 
required  will  be  approximately  two  orders  of  magnitude  less  than  a 
comparable  one-phase  system.  Furthermore,  the  entire  system  is 
essentially  isothermal  which  will  mean  much  smaller  and  lighter 
radiators.  This  system,  however,  has  many  design  problems  that 


must  be  worked  out  before  It  can  be  used.  Therefore,  many  research 
programs  have  boon  started  to  solve  the  problems  associated  with 
two-phase  systems. 

This  study  deals  with  the  design  of  the  evaporator.  More 
specifically,  the  objects  of  this  program  are  to  do  an  in-depth 
review  of  the  forced-convective  boiling  process  to  determine  the 
effect  that  gravity  will  have  on  it,  to  set  up  an  experimental 
evaporator  to  visualize  the  general  performance  and  problems  of 
straigh-tube  evaporators,  to  investigate  methods  of  solving  the 
problems  associated  with  them,  to  review  the  methods  available  for 
calculating  local  beat-transfer  coefficients,  to  develop  a  best 
possible  approach,  to  develop  a  method  for  designing  straight-tube 
evaporators,  and  to  develop  a  computer  program  to  design 
straight-tube  evaporators. 


IT  REVIEW  OF  THE  FORCF.P-CONVECTIVE  BOILING  PROCESS 


In  order  to  understand  tiro  relationship  between  the 
forced-convective  boiling  process  and  gravitational  acceleration,  a 
thorough  study  of  the  boiling  process  was  conducted.  The 
forccd-convecti vc  boiling  process  is  explained  and  the  manner  in 
which  gravitational  acceleration  affects  the  heat-transfer  process 
is  examined  in  this  section, 

GENERAL  DISCUSSION 

The  forced-convective  boiling  process  can  be  broken  down  into 
five  regions:  one-phase  forced  convection;  subcooled  boiling; 
saturated,  nucleate  boiling;  two-phase  forced  convection;  and  the 
liquid-  deficient  region.  This  process  is  shown  in  Figure  1,  and 
the  variation  of  the  heat-transfer  coefficient  associated  with  each 
region  is  shown  in  Figure  2.  The  one-phase,  forced-convection 
region  is  the  initial  condition  normally  encountered  in  an 
evaporator.  It  is  possible,  however,  that  this  region  may  be 
absent.  In  this  region  heat  is  conducted  into  and  is  swept  away  by 
the  fluid. 

As  the  temperature  of  the  wall  increases  along  the  length,  it 
reaches  a  point  where  bubble  micleation  can  begin.  This  is  the 
start  of  the  subcooled-boiling  region.  As  the  wall  temperature 
increases  further,  the  bubbles  grow  larger  and  begin  to  depart  the 
surface.  Since  the  bulk  of  the  fluid  is  still  subcooled,  the 
bubbles  condense  upon  leaving  the  warmer  boundary  layer.  Although 
the  net  amount  of  vapor  generated  in  this  region  is  very  small  ,  the. 
heat-transf er  coefficient  is  very  large. 


■Oft* 


LIQUID-DEFICIENT 

REGION 


TWO-PHASE  FORCED 
CONVECTION  REGION 


SATURATED  NUCLEATE 
BOILING  REGION 


ONE-PHASE  LIQUID 
REGION 


Figure  1:  Regions  associated  with  forced-convective  boiling  in  tubes. 


When  the  bulk  fluid  temperature  reaches  saturation  conditions, 
saturated-nucleate  boiling  occurs.  In  this  region,  normal  nucleate 
boiling  occurs  with  a  corresponding  increase  in  vapor  quality. 

Over  this  region  the  heat-transfer  coefficient  stays  essentially 
constant . 

As  more  vapor  is  produced,  the  liquid  layer  on  the  surface 
becomes  thinner.  Eventually,  this  liquid  layer  becomes  thin  enough 
to  prevent  the  liquid  immediately  adjacent  to  the  wall  from 
superheating  to  the  point  of  allowing  bubble  nucleation.  At  this 
point  the  mechanism  of  heat  transfer  switches  from  boiling  to 
liquid-surface  evaporation.  Since  bubble  nucleation  is  completely 
suppressed  and  the  heat  is  carried  away  by  the  vapor  at  the  liquid/ 
vapor  interface,  this  region  is  called  the  two-phase,  forced- 
convection  region.  This  region  has  the  highest  heat-transfer 
coefficients . 

The  final  region  occurs  when  the  fluid  quality  and  velocity 
increase  to  a  point  where  all  of  the  liquid  becomes  entrained  in 
the  vapor.  At  this  point  the  walls  dry  up  and  the  heat-transfer 
coefficient  decreases  dramatically.  This  "burn  out"  point  must  be 
avoided  at  all  costs  in  evaporator  design.  Heat  transfer  in  the 
liquid-deficient  region  will  not  be  considered  in  this  study. 

In  the  following  text,  two  different  boundary  conditions  will 
be  used  to  further  illustrate  the  boiling  process  in  each  region. 
The  first  boundary  condition  will  be  a  constant  heat  flux  at  the 
evaporator  wall.  The  second  boundary  condition  will  be  a  constant 
evaporator  wall  temperature.  In  reality  a  combination  of  both  will 
be  seen;  but,  for  most  design  problems,  one  of  the  two  boundary 


conditions  can  be  used  with  reasonable  accuracy. 


SUBCOOLED  BOILING 

Subcooled  boiling  occurs  when  the  liquid  adjacent  to  the  wall 
is  superheated  to  a  point  where  bubble  growth  can  occur.  The 

heat-transfer  coefficient  associated  with  this  region  is  much 

s 

larger  than  in  the  one-phase  region.  This  heat-transfer 

coefficient  increases  until  saturated-nucleate  boiling  is  achieved. 

The  heat-transfer  coefficient  begins  to  increase  even  before 

bubble  nucleation  can  be  detected.  A  photographic  study  by  Jiji 
f  31 

and  Clark1  J  demonstrated  this  by  showing  that  a  considerably 
higher  surface  heat  flux  is  required  to  permit  visual  detection  of 
bubbles  than  to  cause  an  improvement  in  the  heat-transfer 
coefficient.  This  Increase  Is  caused  by  hot  Jets  of  fluid  from  the 
boundary  layer  being  projected  into  the  cooler  core  of  the  flow. 

For  the  case  of  a  constant  heat  flux*  the  onset  of  subcooled 
boiling  occurs  in  the  following  manner.  In  the  one-phase,  forced- 
convection  region  the  thermal-boundary  layer  next  to  the  surface 
heats  up  as  more  heat  is  absorbed.  Eventually  a  point  is  reached 
where  bubble  nucleation  can  occur.  A  bubble  will  grow  when  the 
pressure  inside  the  bubble  is  larger  than  the  pressure  outside, 

The  equilibrium  condition  for  the  spherical,  vapor  bubble  of  radius 
r  is: 


where  Py  is  the  pressure  of  the  vapor,  P^  is  the  pressure  of  the 
liquid,  and  o  is  the  surface  tension  of  the  fluid.  The  liquid 
superheat  required  to  maintain  this  equilibrium  can  be  related  to 
the  pressure  difference  with  the  use  of  the  Clausius-Clapeyron 
thermodynamic-equilibrium  relation: 
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(11-2) 


where  AT  is  the  superheat  of  the  liquid  at  the  wall  (T  -T  ). 

S  3t  W  8  At 

On  any  surface  there  are  minor  imperfections  that  contain  a 
trace  of  vapor  or  noncondensable  gas.  These  imperfections  deter¬ 
mine  the  superheat  required  for  bubble  nucleation.  Therefore,  by 
combining  equations  <11—1 )  and  (11-2): 
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it  can  be  seen  that  the  larger  the  size  of  the  imperfection  the 
lower  the  superheat  that  is  required  for  nucleation. 

Hsu1  postualted  that  bubble  nuclei  on  a  heated  wall  will  grow 
only  If  the  lowest  temperature  on  the  bubble  surface  is  greater 


than  the  superheat  required  for  growth.  Earlier  work  by  this 

.(5] 


investigator1 


confirmed  this.  For  a  bubble  of  radius  r£  to  grow 


in  a  uniformly-heated  fluid  of  temperature 


r&  sat 

h.  v 
fg  v 


In  1  + 


(II-4) 


For  a  heated  surface  the  temperature  gradient  can  be  approximated 
in  a  linear  form: ^ 
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(11-5) 


where  T^(y)  is  the  local,  liquid  temperature  and  y  is  the  distance 
away  from  the  wall.  For  a  nucleus  of  radius  r£  to  begin  to  grow 
equations  (11-4)  and  (11-5)  must  be  tangent^.  Therefore: 


dT.(y)  dT. 
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Bergles  and  Rohsenow1  J  obtained  a  graphical  solution  of  equations 
(II-4)  and  (11-5)  subject  to  equations  (II-6)  for  water  over  the 
pressure  range  of  15-2000  psia: 


jlS.S/P0-0231) 


(11-7) 


Davis  and  Anderson  J  developed  an  analytical  solution  of 
equations  (11-4),  (H-5),  and  (11-6): 
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which  works  well  for  water. 

rai 

Frost  end  Dzakovic1  proceeded  in  the  sente  Banner  as  Davis  and 

Anderson;  but*  they  assumed  that  nucleation  occured  when  T^fy) 

2 

equaled  at  a  distance  of  r(Pr)  rather  than  r^.  Thus: 


kS,hf« 


(II-9) 


This  equation  compares  favorably  with  experimental  data  for  a 
variety  of  fluids. 

As  the  temperature  increases  further,  bubbles  form  on  the 
surface,  but  condense  before  they  detach.  During  this  phase  jets 
of  wanner  liquid  continue  to  flow  into  the  cooler  core  by  thermo¬ 
capillarity.  Thermocapillarity  is  caused  by  the  fact  that  the 
surface  tension  varies  around  a  bubble  growing  in  a  temperature 
gradient.  The  resultant  force  pulls  the  interfacial  liquid  film 
outward  (from  hot  to  cold)  which  then  drags  adjacent  warmer  liquid 
with  it.  Cooler  fluid  then  rushes  in  to  replace  the  warmer  fluid. 
Brown1  1  suggested  that  this  phenomenon  can  transfer  a  considerable 
fraction  of  the  heat  flux.  The  turbulence  Induced  by  thermo¬ 
capillarity  explains  why  the  average  temperature  of  the  boundary 
layer  is  below  saturation.  Thermoeapillarlty  also  causes  the 
temperature  of  the  interface  between  the  one-phase  core  and  the 
bubble-laden  boundary  layer  to  fall  with  increasing  heat  flux 
because  of  increased  turbulence. 

In  the  final  stage  of  subcooled  boiling,  when  the  temperature 
approaches  saturation,  the  bubbles  grow  and  detach  from  the 
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surface.  However,  they  still  condense  ns  they  migrate  into  the 
cooler  core  of  the  flow.  The  major  forces  affecting  the  flow 
involved  in  subcooled  boiling  are  the  force  exerted  by  the  fluid 
flowing  through  the  tube  and  the  thermocnpillary  force  at  the 
Chermal-boundary-layer  interface.  Neither  of  these  forces  are 
effected  by  gravitational  acceleration.  (Buoyancy  is  less 
important  than  flow  shear  in  this  regime.)  Furthermore, 
gravitational  acceleration  has  no  effect  on  bubble  nucleation.  For 
these  reasons  the  heat-transfer  coefficient  measured  in  a  1-g 
environment  will  be  the  same  as  that  measured  in  a  0-g  environment 
for  subcooled,  forced-convective  boiling. 

In  actual  practice  the  equations  presented  in  this  section  may 
predict  lower  than  what  is  observed,  especially  at  lower 

pressures  and  low  heat  fluxes.  This  is  due  to  the  fact  that  the 
equations  assumed  a  wide  range  of  activation  sites.  On  metal  or 
glass  tubes  there  may  not  be  large  enough  nucleation  sites  for 
boiling  to  occur  even  though  the  equations  predict  it. 

For  the  case  where  there  is  a  constant  wall  temperature, 
subcooled  boiling  will  begin  immediately.  However,  the  boiling 
process,  discussed  above,  will  be  the  same. 

SATURATED-NUCLEATE  BOILING 

When  the  bulk-fluid  temperature  reaches  saturation  conditions, 
the  vapor  bubbles  no  longer  condense.  Instead  they  continue  to 
grow  and  coalesce  forming  even  larger  bubbles.  Therefore,  the 
vapor  quality  begins  to  increase. 

In  this  region  the  heat-transfer  mechanism  is  very  similar  to 
nucleate-pool  boiling.  In  fact,  the  flow  rate  of  the  fluid  has 


very  little  influence  on  the  heat-transfer  coefficient^. 

Assuming  that  the  heat  transfer  is  caused  by  two  additive  effects: 
forced-convection  and  nucleate  boiling;  the  forced-convectien 
component  can  be  Ignored.  Therefore,  an  equation  for  the 

heat-transfer  coefficient  in  nucleate,  pool  boiling  can  be  used  to 

> 

calculate  the  heat-transfer  coefficient  for  forced-convection, 
saturated,  nucleate  boiling.  This  has  been  done  successfully  by 
several  investigators^^'  However,  one  must  be  careful  in 

this  approach  because  subcooling  greatly  affects  pool  boiling, 
while  having  some  subcooling  (the  average  temperature  is  the 
saturation  temperature)  in  the  liquid  core  for  forced  convection 
has  no  effect. 

The  reasons  for  the  large  increase  in  the  heat-transfer  coeffi¬ 
cients  for  nucleate  boiling  have  been  debated  for  years  with  no 
concise  agreement The  main  mechanism,  however,  seems  to  be 
the  turbulence  caused  by  bubble  generation  and  departure.  Before 
the  bubble  grows,  a  thermal-boundary  layer  is  built  up.  When  the 
temperature  near  the  surface  gets  high  enough  for  bubble  growth, 
the  bubble  grows  outward,  pushing  the  boundary  layer  with  it. 
Furthermore,  when  the  bubble  detaches  from  the  surface  it  then 
carries  this  thermal-boundary  layer  into  the  bulk  of  the  fluid, 
thus,  directly  transporting  energy  outward.  Cooler  fluid  then 
rushes  back  to  the  surface  to  complete  the  cycle. 

The  major  factors  affecting  heat  transfer  during  saturated- 
nucleate  boiling^  are  the  heat  flux  and  the  pressure.  The  flow 
rate,  quality,  and  core  subcooling  have  little  effect  in  this 
region.  Furthermore  the  heat-transfer  coefficient  remains  constant 
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throughout  the  region  because  of  the  constant  bulk  temperature, 
(equal  to  the  saturation  temperature).  And,  the  boundary 
conditions  at  the  wall  have  no  effect. 

In  this  region  the  major  forces  acting  on  the  fluid ^  are 
surface  tension,  inertia,  bouyancy,  drag,  and  the  Bernoulli  effect. 
Surface  tension  holds  the  bubble  to  the  wall  while  the  inertia  of 
the  liquid  surrounding  a  growing  bubble  attempts  to  pull  it  from 
the  wall.  The  drag  force  caused  by  the  fluid  flowing  around  the 
bubble  attempts  to  shear  it  from  the  wall.  The  bouyancy  force 
pulls  the  bubble  upward,  pulling  it  from  the  wall  or  pushes  it  onto 
the  wall  depending  on  the  surface  orientation.  Upon  departure  from 
the  surface,  bouyancy  makes  the  bubble  rise  while  the  Bernoulli 
effect  created  by  the  velocity  gradient  within  the  tube  causes  the 
hubbies  to  migrate  to  the  center.  As  long  as  the  flow  is  not 
stratified,  the  differences  between  the  heat-transfer  coefficient 
in  1-g  and  0-g  for  a  round  tube  should  be  small.  This  is  because, 
in  l-g»  gravity  helps  on  the  bottom  of  the  tube  by  helping  to 
remove  the  vapor  but  hinders  on  the  top  by  keeping  the  vapor  near 
the  surface.  The  average  effect  should,  therefore,  be  similar  to 
0-g  conditions. 
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TWO-PHASE,  FORCED  CONVECTION 

As  the  vapor  quality  increases  the  vapor  migrates  towards  the 
center  creating,  at  first,  a  slug-flow  pattern  which  is  followed  by 
an  annular-flow  pattern.  A  slug-flow  pattern  is  where  large 
bubbles  or  "slugs"  of  vapor  are  separated  from  one  another  by  thin. 
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liquid  regions.  As  the  vapor  quality  increases  the  liquid  regions, 

separating  the  slugs,  disappear  leaving  an  annular-flow  pattern 

where  the  vapor  core  is  separated  from  the  wall  by  a  thin  liquid 

layer.  Flow  patterns  and  the  conditions  for  their  appearance  are 

discussed  extensively  in  Reference^. 

Eventually,  a  point  is  reached  where  the  liquid  layer  adjacent 

to  the  wall  is  thin  enough  that  all  of  the  heat  flux  is  conducted 

through  the  liquid  without  nucleate  boiling.  The  mechanism  of  heat 

transfer  then  switches  from  nucleate  boiling  to  evaporation  at  the 

liquid/vapor  interface.  This  regime,  therefore,  is  called  the 

two-phase,  forced-convection  region. 

In  order  to  suppress  bubble  nucleatlon,  the  effective  thermal 

conductivity  of  the  liquid  layer  must  be  high  enough  so  that  the 

liquid  cannot  be  superheated  to  a  point  where  bubble  growth  occurs. 

The  equations  used  earlier  to  predict  the  onset  of  subcooled, 

nucleate  boiling  are  applicable  here  too.  As  mentioned  earlier, 
f  8 1 

Frost  and  Dzakowic1  J  predicted  the  bubble  growth  would  begin  when: 
k.h  /AT  \2 

^ONB  *  6ot  v7  V  Pr  y  (II -9) 

sat  fg>  “  ' 

Combining  this  equation  with  Newton's  law  of  cooling: 

4SNB  *  h0TPFC^Tsat  (11-10) 


where  is  the  heat-transfer  coefficient  for  the  onset  of 

two-phase,  forced  convection,  and  qg^g  is  the  heat  flux  for  the 


suppression  of  nucleate  boiling.  Setting  q 
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qgNB  an  expression 
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for  h.-*,.,.,  for  the  onset  of  two-phase,  forced  convec- 
UTr  FC 

tlon  can  be  found. 


k£*VTsat 

OTPFC  “  *  w^2 


(11-11) 


This  equation  for  is  t*ie  point  where  bubble  nucleation 

is  suppressed.  As  in  subcooled  boiling,  this  number  may  be  lower 
than  what  is  actually  found  in  practice  due  to  the  absence  of  large 
nucleation  sites. 

In  this  two-phase,  forced-convection  region  the  heat  is  trans¬ 
ferred  across  the  liquid  layer  and  evaporates  liquid  at  the  liquid/ 
vapor  interface.  Therefore,  the  heat-transfer  coefficient  is  very 
dependent  on  the  thickness  of  the  liquid  layer.  This  thickness  is 
governed  by  the  quality.  Therefore,  h^p^  is  very  dependent  on 
quality.  Furthermore,  as  the  thickness  of  the  liquid  layer 
decreases,  the  temperature  drop  across  the  liquid  layer  decreases. 
Because  of  this^  extremely  high  heat-transfer  coefficients  can  be 
achieved. 

Another  factor  that  affects  the  heat-transfer  coefficient  is 
the  mass  flux^.  Higher  mass  fluxes  increase  the  vapor  velocity 
which  increases  the  evaporation  at  the  liquld/vapor  interface. 
Unlike  the  saturated-nucleate-boiling  regime,  the  heat  flux  has  no 
effect  on  the  heat-transfer  coefficient. 

Since  the  heat-transfer  coefficient  increases  along  the  tube 
with  increasing  quality,  the  wall  temperature  of  the  tube  will 
decrease  along  the  tube  for  the  case  of  a  constant  heat  flux.  For 
the  case  of  a  constant  wall  temperature,  the  heat  flux  will 
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increase  along  the  tube. 

In  this  regime,  the  only  contribution  of  gravity  to  the  heat 

transfer  process  is  in  the  distribution  of  film  thickness  around 

the  circumference  of  the  tube.  However,  the  average  value  will  be 

[51 

the  same  ns  the  constant  value  found  under  0-g  conditions1  ,  if 
the  fluid  wets  the  tube.  Therefore,  heat  transfer  in  this  regime 
should  be  relatively  insensitive  to  gravity.  Furthermore,  all 
analytical  treatments  of  this  regime  have  used  a  constant,  fluid 

thickness  around  the  perimeter  (0-g  condition)  successfully  in  past 

,  [5] 
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LIQUID  DEFICIENT  REGION 

As  the  quality  increases  in  the  two-phase,  forced-convection 
region  the  liquid  layer  getB  thinner,  while  the  vapor  velocity 
increases.  Eventually  a  point  is  reached  where  the  liquid  becomes 
totally  entrained  in  the  vapor  core,  thus  causing  the  walls  to 
dryout.  When  dryout  conditions  are  reached,  the  heat-transfer 
coefficient  decreases  drastically.  With  n  constant  heat  flux,  this 
causes  the  wall  temperature  to  increase  tremendously.  Often  it 
will  increase  above  the  melting  point  of  the  evaporator  walls,  thus 
"burning  out"  the  evaporator.  For  the  case  of  a  constant  wall 
temperature,  the  heat  flux  decreases  drastically.  This  causes  the 
liquid  deficient  region  of  the  evaporator  to  be  essentially 
useless.  This  condition,  therefore,  must  be  avoided. 
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Ill  FLOW-VISUALIZATION  EXPERIMENT 
In  order  to  observe  the  general  operation  of  and  problems 
associated  with  straight-tube  evaporators,  a  flow-visualization 
experiment  was  set  up.  The  experimental  apparatus  and  the  results 
of  the  experiment,  are  discussed  below. 


EXPERIMENTAL  APPARATUS 

The  flow-visualization  loop  consisted  of:  an  evaporator,  a 
coridenser/accumulator,  a  preheater,  a  pump,  a  flow  meter,  and  an 
adjustable  orifice.  This  set  up  is  shown  in  Figure  3.  The 
evaporator  consisted  of  a  long,  small -diameter,  glass  tube  and  will 
be  discussed  in  greater  detail  later.  The  condenser/accumulator 
consisted  of  a  copper  coil  (condenser)  submerged  in  a  bucket  of 
water  (accumlator).  The  steam  and  water  mixture  from  the  evapora¬ 
tor  flowed  through  tha  copper  coil  and  emptied  into  the  bottom  of 
the  bucket.  The  copper  coil  thus  served  as  a  condenser  and  as  a 
heater  for  the  water  supply  in  the  accumulator.  The  accumlator  was 
open  to  the  atmosphere.  The  preheater  was  a  hot  plate  beneath  the 
accumulator.  It  was  used  to  help  raise  the  temperature  of  the 
water  to  near  saturation.  The  flow  meter,  a  rotameter,  consisted 
of  a  small,  verticle,  glass  tube.  A  steel  and  a  rubber  ball  rose 
within  the  tube  with  an  increase  in  the  flow  rate.  The  flow  meter 
also  contained  a  valve  that  could  be  used  to  regulate  the  flow. 

The  purpose  of  the  adjustable  orifice  will  be  discussed  later. 


Instrumentation 

Total  instrumentation  consisted  of  three  pressure  meters, 
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three  thermocouples ,  an  ammeter,  and  a  voltmeter.  Two  of  the 
pressure  meters  were  set  up  to  measure  the  pressure  drop  across  the 
evaporator  and  across  the  adjustable  orifice.  The  third  pressure 
gage  measured  the  pressure  of  the  flow  prior  to  entering  the 
orifice.  The  three  thermocouples  were  located  before  and  after  the 
evaporator  and  in  the  accumulator.  The  ammeter  and  voltmeter  were 
used  to  measure  the  power  input  to  the  evaporator. 


Evaporator 

The  evaporator  consisted  of  a  long,  thin,  glass  tube.  Heat  was 
supplied  to  the  tube  by  nichrome  heating  wire  which  was  coiled 
around  the  tube.  The  heat  from  the  wire  was  conducted  through  the 
glass  tube  which  caused  the  water  to  boil  as  it  passed  through  the 
tube. 

The  purpose  of  experiment  was  to  visualize  the  boiling  process 
in  a  spacecraft  evaporator.  To  do  this  it  was  necessary  to 
minimize  the  influence  of  gravity.  On  earth,  gravity  tends  to 
dominate  two-phase  flow  through  tubes,  while  in  space,  surface 
tension  tends  to  dominate.  However,  in  very  small-diameter  tubes, 
surface  tension  dominates  over  gravity.  For  thi6  experiment,  this 
point  was  assumed  to  be  the  point  where  the  liquid  is  able  to  wick 
itself  completely  around  the  tube  so  that  a  stratif ied-flov  pattern 
cannot  exist.  From  earlier  work^  this  point  is: 
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Using  water  properties: 
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f  2(5, 89X10 ‘"‘NT/b)  | 

r  “[  (955kg/R3)(9.8si/s2)J 

t  *  0.0035&  -  O.llft  »  l/8in 

Therefore,  tubes  with  *s  inch  diameters  or  less  were  used  for  the 
evaporator. 


RESULTS  OF  FLOW  VISUALIZATION 

The  flow  entered  the  evaporator  subcooled.  Almost  immediately 
very  small  bubbles  would  appear  at  the  evaporator  wall  (subcooled 
boiling).  Although  these  bubbles  remained  small,  those  generated 
further  along  the  tube  did  grow  as  the  bulk  temperature  increased. 
They  also  departed  from  the  surface  and  recondensed  in  the  flow. 

A3  the  bulk  temperature  of  the  flow  approached  saturation,  the 
bubbles  began  to  grow  larger  and  coalesce.  When  the  bulk 
temperature  reached  saturation  the  flow  pattern  changed  from  bubbly 
flow  to  slug  flow.  In  this  slug-flow  regime  the  liquid  layer 
between  the  vapor  core  and  the  wall  was  thin  enough  to  suppress 
nucleate  boiling.  Still  further  down  the  tube,  the  slugs  expanded 
into  an  annular-flow  pattern.  This  pattern  established  in  the 
evaporator  was  always  the  same  regardless  of  power  input  or  flow 
rate.  The  only  difference  being  whether  enough  power  was  being 
supplied  to  reach  the  flow  quality  necessary  for  annular  flow. 

Although  the  flow  pattern  established  in  the  evaporator  was 
consistent,  it  was  not  stable.  This  instability  was  caused  by  the 
extreme  pressure  oscillations  associated  with  slug  flow.  In  a 
spacecraft  evaporator  these  pressure  oscillations  cannot  be 
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allowed.  They  would  cause  vibration,  disruption  of  flow 
uniformity  in  paratJ el-flow  networks,  wear  on  pumps,  and 
instabilities  in  the  condenser.  Therefore,  before  straight-tube 
evaporators  can  be  used,  a  means  of  stabilizing  the  flow  must  be 
devised. 

In  an  effort  to  solve  this  problem,  an  orifice  was  installed 
immediately  in  front  of  the  evaporator  tube.  If  the  orifice  could 
be  sized  so  that  the  pressure  drop  through  the  evaporator  is  much 
smaller  than  the  pressure  drop  through  the  orifice,  the  flow  should 
be  stable  in  front  of  the  evaporator.  This  would  eliminate  the 
problems  of  parallel-flow  instability  and  pump  wear.  To  test  this 
theory,  an  adjustable  orifice  was  placed  ahead  of  the  evaporator. 
The  pressure  was  measured  at  point  1  (see  Figure  3) ,  and  the 
pressure  drop  was  measured  between  points  1  and  2  and  between  2  and 
3.  The  results  of  this  test  were  very  conclusive  -  this  method 
will  not  work!  The  orifice  in  this  experiment  was  varied  until  the 
pressure  drop  across  the  orifice  was  approximately  an  order  of 
magnitude  greater  than  the  pressure  drop  through  the  evaporator. 

And  still,  the  flow  oscillations  could  not  be  damped  out  prior  to 
the  orifice.  The  pressure  readings  are  shown  in  Table  I  for 
various  mass  flow  rates.  An  investigation  of  the  accumulator 
showed  the  flow  to  be  steady  to  the  pump ,  therefore  the  only  cause 
of  these  instabilities  could  be  the  slug  flow.  Any  further 
increase  in  the  pressure  drop  across  the  orifice  was  deemed 
unacceptable  because  the  corresponding  Increase  in  pump  power 
required  would  eliminate  the  major  advantage  of  the  two-phase 
systems  having  a  low  pump-power  requirement . 
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TV  TWO-PHASE,  HEAT-TRANSFER  EQUATIONS 

In  order  to  determine  a  method  of  calculating  the  heat-transfer 
coefficients  for  forced-convective  boiling,  a  review  of  past  work 
is  presented  here.  The  information  obtained  in  this  section  ia 
then  used  to  develop  a  better  method  in  Section  V. 

Many  empirical  equations  have  been  developed  to  predict  the 
heat-transfer  coefficients  in  forced-convection-boiling  heat 
transfer.  Most  of  the  equations  were  developed  for  only  one  region 
of  boiling,  while  others  use  weight  factors  with  several  equations 
to  make  a  general  equation.  The  weight  given  to  a  certain  equation 
depends  on  the  local  flow  conditions.  This  section  describes  some 
of  the  newer  and/or  more  popular  equations. 

One  of  the  easiest  equations  to  use  is  that  of  McAdams  et 
r  121 

al.  They  developed  an  empirical  equation  for  the  heat  flux  for 

water  during  subcooled,  nucleate  boiling: 

9  »  C'AT  3-86  (1V-1) 

fiat 

where  C'  is  a  dimensionless  constant  that  ranges  from  0.190  to 
0.074  depending  on  the  extent  of  degassing  of  the  water  and  T 
is  the  difference  between  the  wall  temperature  and  the  saturation 
temperature.  They  found  q  to  be  Independent  of  water  velocity, 
pressure,  degree  of  subcooling,  and  diameter,  The  equation  is  good 
for  pressures  of  30  psia  to  90  psia.  The  main  problem  with  this 
equation  is  that  it  was  developed  solely  for  water  and  is  not 
applicable  to  other  fluids. 


As  mentioned  previously,  pool-boiling  equations  have  been  used 


successfully  for  calculating  the  heat-transfer  coefficients  in  the 


subcooled  and  saturated-nucleate-boiling  regions.  Specifically, 

Lnvin  and  Young^^,  and  Myers  and  ICatsJ*^  have  successfully  used 
r  131 

McNelly1 s  equation1  1  for  the  nucleate,  pool-boiling  heat-transfer 
coefficient  to  correlate  their  saturated,  forced-convection 
nucleate-boiling  data.  McNelly* s  equation  is: 


(1V-2) 


where  C  was  given  as  .225  by  McNelly  but  modified  to  .25  for  the 
inside  of  a  plain  tube  by  Lavin  and  Young.  And: 
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McNelly  developed  this  equation  by  first  looking  very  closely 
at  the  boiling  process,  and  then  doing  a  dimensional  analysis.  He 
also  used  the  heat  flux  as  the  driving  mechanism  rather  than  the 
temperature  difference.  It  should  be  noted  that  like  the  McAdams 
equation,  the  diameter  has  no  effect  (divides  out)  and  the  mass 
flux  has  no  effect  (developed  for  pool  boiling).  This  equation  is 
good  for  all  fluids. 

Lavin  and  Young  also  went  on  to  develop  an  equation  for  the 
two-phase,  forced-convection  region.  To  correlate  their  data  they 
developed  the  equation: 


(IV-4) 
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where  h^,  is  the  two-phi'  forced-convection,  heat-transfer 
coefficient;  h^p  is  the  one-phase,  heat-transfer  coefficient  for 
the  liquid  phase;  G  is  the  mass  flux;  and  x  is  the  quality.  For 
h^p  they  used  the  Sieder-Tate  equation: 

bifB0'023iRei0'8pr°‘33^J  (iv- 5 

where  the  Reynolds  number  is  defined  as: 


_  Gd( 1-X 

%  y  (IV-6) 


Since  this  equation  was  developed  for  the  two-phase,  forced- 
convection  region,  it  is  very  dependent  on  the  quality  and  relative¬ 
ly  independent  of  the  heat  flux. 

Lockhart  and  Martinelli  devised  a  parameter,  X,  to  correlate 
the  pressure  drop  for  two-phaBe  flow  in  tubes.  This  parameter  is 
also  used  to  calculate  the  liquid-film  thickness  in  two-phase 
flows.  Since  the  heat-transfer  coefficient  is  very  dependent  upon 
the  liquid  layer  adjacent  to  the  surface  in  two-phase,  forced 
convectionj many  investigators  have  used  X  to  successfully  correlate 
their  data.  They  used  an  equation  of  the  form: 

(IV-7) 
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where  Is  the  Loekhart-Martinelli  parameter  for  turbulent 
liquid,  turbulent  vapor: 


(IV-8) 


Dengler  and  Addome 


[1M 


found  that: 


(IV-9) 


where  h^  Is  the  heat-transfer  coefficient  for  liquid  alone  found 
by  setting  x-0.0  and  using  equation  (IV-5).  Guerrierl  and 
Talty^"^  found  thats 


(1V-10) 


where  hTT>  is  the  one-phase,  heat-tranafer  coefficient.  Bennett  et 
Lr 

ait  found  that: 


(XV-11) 
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All  of  the  above  correlations  were  compared  by  Chen^^.  He 
found  that  Dengler  and  Addons'  equation  predicted  values  too  large 
when  compared  vlth  other  experimenters  data.  Guerrlerl  and  Xalty's 
equation  predicted  lover  coefficients  than  Dengler  and  Addons'  but 
still  lacked  accuracy,  and  Sennet’s  predicts  hTp  veil  for  vater 
but  not  for  other  fluids.  For  these  reasons  Chen  developed  his  own 
method  which  covers  both  the  saturated,  nucleate-boiling  region  and 
the  two-phase,  forced-convection  region.  He  assumed  that  the  two 
mechanisms,  nucleate  boiling  and  forced  convection,  are  present  to 
6ome  extent  in  both  regions.  Therefore,  he  devised  a  correlation 
using  superposition: 


hTP  “  hNB  +  hFC  (IV-12) 

Where  is  the  contribution  due  to  nucleate  boiling  and  hpC  is 
the  contribution  due  to  forced  convection.  A  Dittus-Boelter  type 
equation  was  used  for  hpr. 


0.023Re 
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(IV-13) 


Chen  argued  that  k^p  and  PrTp  can  be  equated  with  the  .liquid 
values.  Furthermore,  he  defined  a  parameter  F  as: 


f  » 


(IV-H) 
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The  only  unknown  In  hpC  then  is  F.  However,  since  F  Is  a  flow 
parameter  it  can  be  related  to  Xt{;.  Chen  used  Foster  end  Zubers' 
analysis  of  nucleate  boiling  for  h^g. 
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He  also  defined  a  nucleate-boiling,  suppression  factor,  S,  which 
relates  the  mean  superheat,  AT,  to  the  wall  superheat,  At  j 

8®C 


S 


(1V-16) 


using  the  Clausius-Clapeyron  equation: 
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Combining  equations:  (XV— 1 5 ) ,  (IV— 16)  and  < IV— 17);  h^g  equals: 


*  0.00122 


v  0.79-  0. 1*5  0.1*9  0.25.  0.1*9 

k*  "a..  _ 


0.5,.  0.29n  0.21*_ 

0  U£  pv  Tsat 


0.75  0.75 

Tfg 


ATS  ( IV— 18) 


h  “  hN3  +  ^FC 


(IV“19) 


The  Faraaeters  S  and  F  are  graphed  in  Figures  4  end  5.  The  Chen 
equation  has  been  correlated  with  the  data  of  mary  investigators 
with  85%  of  the  data  being  within  -  15. IX.  The  problem  with  the 
Chen  equation  is  that  it  is  not  easy  to  program  on  a  computer  due 
to  the  use  of  charts  to  find  S  and  F. 

A  final  method  for  calculating  Hjp  was  proposed  recently  by 

fig] 

BJorge,  Hall,  and  Rohsenow.1  J  They  also  used  a  superposition 
formula  where: 


1Tp  *  1FC  +  qNB  -  qIB 


(1V-20) 


for  the  high  quality  region  (x  >  0.05),  and 
qTP  “  [qFC  +  ^  qN3  “  qlB^  ] 1 


(1V-21) 


for  the  low  quality  region.  qfB  is  the  heat  flux  at  the  incipient 
boiling  point.  Since  q^g  and  q^  are  proportional  to  T^: 


^TP  =  4FC  + 


\  J  . 


(1V-22) 


for  the  high  quality  region,  and  for  the  low  quality  region: 
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Figure  4;  Reynolds  number  factor  F  versus  the  Lockhart- 
Martlnelli  parameter  Xtt. 
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Figure  5: 


Suppression  factor  S  versus  the  Reynolds  number. 
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For  the  high  quality  region ,  the  equation  developed  by  Traviss  et 
al^*^.  and  modified  by  Hall^^  was  used  for  q^: 

*FC  “  - - “O - - - “Aisat  (IV-24) 


where : 
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F2  »  5Pr£+51n(l+5?r£)+2.51n(0.0031Re£0,8:i2)  (Re>1125)  (1V-26*) 
F?  -  5Pr£+51n[l+Pr£(0.096^Re£D,585-l)]  (50<Re<1125)  (IV-26b) 

■  0.0707Pr£Re£0,5  (Re>50)  <IV"26c) 
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The  equation  was  developed  by  Ttaviss  et  el.  as  a  correlation  for 
condensation  in  annular  flow.  Since  two-phase,  forced-convection 
evaporation  is  the  opposite  of  forced-convection  condensation,  this 
equation  should  work  for  two-phase  evaporation.  Hall  found  this  to 
be  true,  with  a  slight  aodification  of  the  coefficient. 

For  the  subcooled  and  low-quality  region,  Colburn's  equation 
was  uaed: 


(1V-27) 


where  the  subscript  f  means  that  the  property  should  be  evaluated 
at  the  film  temperature,  (Tw  +  /2 »  the  eubscript  b  refers  to  the 

bulk  temperature,  end  qF{.  can  be  found  with  equation  (IV-3) 

The  nucleate  boiling  heat  flux  for  both  regions  is  calculated 
by  using  the  Mikic-Kohsenov  correlation: 
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where  is  1.89X10  in  SI  units  for  forced-convection  boiling  of 
water. 

For  the  wall  superheat  at  the  incipience  of  boiling,  they  used 


the  equation  derived  by  Davis  and  Anderson  for  the  high  quality 
region: 


In  the  subcooled  and  lev  quality  region,  this  equation  is  modified 
to  account  for  the  fact  that  there  may  not  be  large  nucleation 

sites : 


ATsat,IB  *  l-n(i*rN  "  K*Tsc)  (rtang>rmax^  <IV”30> 

where  rtang  *s  the  radius  of  tangency  between  equation  11-4  and 
11-5  in  the  subcooled-boiling  section  or: 
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Furthermore : 
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Where  t  *  10  Am. 
max 

It  should  be  noted  that,  in  the  method  devised  by  Bjorge,  Mall 
and  Rohsenov,  hy^  for  the  high-quality  region  is  completely 
different  from  hy^,  in  the  low-quality  and  subcooled  regions.  In 
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the  low-quality  and  subcooled  region  represents  the  one-phase, 
forced-convection,  heat-transfer  coefficient,  equation  (XV-27),' 
while  in  the  high-quality  region,  hpc  represents  the  heat  transfer 
coefficient  for  forced  convection  with  a  phase  change  at  the 
boundary,  equation  (IV-24) . 


V  PROPOSED  METHOD  FOR  DESIGN  OF  STRAIGHT-TUBE  EVAPORATORS 

In  this  section  the  information  that  was  presented  in  the 
proceeding  three  sections  Is  combined  to  develop  a  new  approach  to 
designing  straight-tube  evaporators,  and  a  new  method  for  sizing  the 
length  of  the  evaporator  tube  is  developed.  Next,  a  new  way  of 
calculating  the  local  heat-transfer  coefficients  and  how  to 
stabilize  the  flow  is  presented.  This  is  followed  by  a  way  to 
calculate  the  pressure  drop.  And  finally  all  of  this  information 
was  used  to  develop  a  computer  design  program  which  is  also 
explained  in  this  section. 

DESIGN  METHOD 

Presently,  the  length  of  the  tube  necessary  for  a  certain  size 
of  evaporator  is  calculated  by  treating  the  one-phase  and  two-phase 
sections  separately.  First  the  length  of  the  one-phase  section  is 
calculated  by  using  the  one-phase  heat-transfer  coefficient,  the 
energy  it  takes  to  reach  saturation  conditions,  and  the  desired 
wall  temperature.  This,  of  course,  gives  a  certain  heat  flux  in 
the  one-phase  section. 

After  the  calculations  are  finished  for  the  one-phase  section, 
the  length  of  the  two-phase  section  is  calculated  in  the  same  way, 
using  appropriate  two-phase  equations.  For  this  calculation  it  is 
assumed  that  boiling  begins  when  the  bulk  temperature  reaches 
saturation.  The  total  length  can  then  be  calculated  by  adding  the 
lengths  and  multiplying  by  a  safety  factor. 

The  problem  with  this  approach  is  that  by  breaking  up  the  tube 
into  two  separate  parts,  the  heat  flux  is  considerably  different  in 


the  two  sections.  In  fact,  since  Q/A  is  proportional  to  h  and  hQp 
is  less  than  h^p,  Q/AQp  will  be  much  less  than  Q/A^.  This, 
however,  will  not  be  the  case  in  an  actual  evaporator.  Instead, 
there  will  nearly  be  a  constant  heat  flux,  such  as  with  electronics 
cooling,  or  a  constant  wall  temperature,  such  as  with  exchanging 
heat  with  a  heat  pipe.  Furthermore,  with  a  constant-  wall 
temperature,  subcooled  boiling  will  occur  rather  than  one-phase 
heat  transfer.  With  a  constant  heat  fluxt the  wall  temperature  will 
vary  with  the  heat-transfer  coefficient,  h,  while  with  a  constant 
wall  temperature,  the  local  heat  flux  will  vary  with  h. 

Because  of  the  inaccuracies  of  the  past  approach,  the  following 

two  methods  of  "sizing"  an  evaporator  have  been  developed.  These 

methods  are  depicted  in  Figures  6  and  7  for  a  constant  heat  flux, 

and  for  a  constant  wall  temperature,  respectively.  For  a  constant 

heat  flux,  first  calculate  the  incipient-boiling  point.  This  is 

the  point  where  the  wall  temperature  is  high  enough  above  the 

saturation  temperature  that  boiling  can  begin.  Next,  calculate  h^p 

and  hTp,  hTp  will  vary  along  the  length  of  the  two-phase  section. 

Using  the  allowable  At  -  T  -T  ,  calculate  the  length  using 

Qtotai*  avera8e  h-pp,  d  and  T  .  This  is  the  length  of  tube 

necessary  for  the  evaporator.  The  length  of  the  subcooled  section 

is  the  length  necessary  to  get  AT  j-  ATg(,  ■  ATJB,  where  AT  j 

■  T  -  T,  ,  AT  ■  T  _  -  T,  ,  and  aTt_  ■  Trn  -  T  From  the 

w  b  sc  sat  b  IB  IB  sat 

incipient  boiling  point  to  the  saturation  point,  subcooled  boiling 
will  take  place  and  h^p  must  be  used.  Prior  to  that  point  h^p  can 
be  used  to  determine  the  wall  temperature. 
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For  a  constant  wall  temperature,  the  heat  flux  will  vary  along 
the  tube.  Since  the  wall  temperature  will  be  higher  than  the 
Incipient-boiling  point,  boiling  will  start  immediately.  There¬ 
fore,  using  appropriate  two-phase  equations,  hTp  and  the  heat  flux 
can  be  calculated.  And  finally,  the  length  can  be  found. 


CALCULATION  OF  THE  HEAT-TRANSFER  COEFFICIENT; 

To  calculate  the  heat-transfer  coefficient,  the  superposition 
method  of  Rohsenow^-18^  was  employed.  This  method  was  used  because 
of  its  success  In  the  past  and  its  adaptability  to  a  computer 
solution.  This  method  assumes  that  the  forced-convection  boiling 
heat  flux  can  be  found  by  summing  the  heat  fluxes  for  pool  boiling 
and  forced-convection.  More  specifically,  for  subcooled  boiling 
flow  ( x=  0 ) : 


(IV-23) 


and  for  saturated-boiling  flow  (x  >  0): 


^TF  “  qFC  +  qNB 
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and  for  high-quality  flows  where  nucleate  boiling  is  suppressed 


^TP  s  qFC 


(V-l) 
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Figure  7;  Flow  chart  for  calculating  the  length  of  an  evaporator 
with  a  constant  wall  temperature. 
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Where  q^  for  the  subcooled  regions  is  the  one-phase, 
forced-convection  heat  flux,  equations  (IV-27)  and  (V-2),  and  qFC 
Is  the  forced-convection,  evaporative  heat  flux,  equation  (IV-24), 
for  the  saturated  regions.  The  force-convective  heat  flux  for  the 
subcooled  region  can  be  calculated  by  Colburn's  equation: 


Kn  •  0 

FC 


0.333 


(IV-27) 


with: 


qFC  *  hFC(ATB«t  4  *rBC) 


(v-2) 


qpc  for  the  high-quality  region  can  be  found  by  Travlss' 
equation: 

qF c  F0d  AT8ftt 


( 1V-24) 
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F2  *  5PV51n(l+5Prji)+2-51n(0.0031Pe£0'812)  (Fe>ll25)  OV-26a) 
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(50<Re<1125)  (IV-26b) 
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for  both  equations  can  be  found  by  using  McNelly's  equation: 


0.31,^0.31 
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where  aTp  is  the  total  heat  flux  and: 


ins  ®  8at  I 
And  finally,  AT$at  can  be  found  with  Frost  and  Dzakowic's 
equation: 


(V-3) 
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where  hpc  is  found  from  equation  (IV-27). 


COMPARISON  WITH  DATA 

The  method  outlined  above  was  compared  with  data  from  Cheng, 
et.alJ2^;  Latsh,  et.alJ^  and  Wright^^.  Cheng's  data  was  for 
subcooled  boiling  water  and  is  shown  in  Figure  8.  In  this  figure 
both  the  method  outlined  above  and  the  original  method  by  Bjorge, 
Hall  and  Rohsenow  is  shown.  Notice  that  there  is  very  little 


boiling.  Data  front  Reference  21. 


difference  between  the  two.  This  is  due  to  the  fact  that  both 
methods  are  the  same  except  for  the  equation  used  to  calculate 
qNg.  For  water  both  these  equations  give  essentially  the  same 
answers.  However,  the  equation  used  by  Bjorge,  Hall,  and  Rohsenow 
is  only  applicable  to  water,  where  the  new  method  developed  in  this 
study  is  applicable  to  any  fluid. 

Latsh's  data  for  one-phase,  forced-convection  and  subcooled 
boiling  is  shown  in  Figure  9.  This  figure  shows  where  the 
calculated  incipient  boiling  point  is.  Notice  that  this  point  is 
where  the  data  makes  a  bend  upward.  This  upward  bend  is  also 
evident  in  the  computer  generated  lines. 

Wright's  data  for  low-quality  and  high-quality  flow  is  shown  in 
Figures  10  through  13.  The  accuracy  of  the  low-quality  predictions 
clearly  shows  that  convective  evaporation  is  taking  place  under 
low-quality  conditions.  In  fact,  in  Appendix  A  where  a  sample 
calculation  is  shown,  forced-convective  evaporation  comprises  all 
of  the  heat  transfer  at  qualities  lower  than  five  percent  under  the 
flow  condition. 

It  should  be  noted  that  this  method  is  similar  to  that  develop¬ 
ed  by  Bjorge,  et.  al .  They  chose  X  =  0.05  as  the  boundary  between 
equation  (IV-23)  and  (IV-22)  and  did  not  use  equation  (V-l). 
However,  the  flow  visualation  experiment  done  in  this  study  and  a 
comparision  of  the  equations  with  Wright's^^  data  showed  that  the 
boundary  between  equation  (IV-23)  and  (IV-22)  should  be  between 
zero  and  one  percent  quality.  For  computational  ease,  this  value 
was  chosen  as  zero  percent  quality.  Furthermore,  sample 
calculations  (see  Appendix  A)  showed  that  nucleate  boiling  Is 


suppressed  at  relatively  low  qualities.  If  equation  (IV-22) 
continues  to  be  used  for  suppressed,  nucleate-boiling  conditions 
the  nucleate-boiling  tens  is  actually  being  subtracted  off.  This, 
of  course.  Is  inaccurate.  Therefore  equation  (V-l)  must  be  used 
under  suppressed,  nucleate-boiling  conditions. 

FLOW  STABILIZATION 

A  method  of  stabilizing  the  flow  is  to  prevent  slug  flow  from 
occurring.  This  can  be  done  by  inducing  an  arcifical,  body  force 
so  that  stratified  flow  occurs  instead  of  slug  flow.  The  simplest 
way  of  accomplishing  this  Is  to  induce  swirl  in  the  flow  so  that 
there  is  a  large  centrifugal  force.  To  Induce  swirl  a  twisted  tape 
Inside  the  tube  can  be  used.  If  the  vicking  height  of  the  fluid  is 
used  as  a  transition  point  between  slug  and  stratified  flow^: 
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The  tangential  acceleration  due  to  fluid  swirl  is: 


(V-6) 
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flux,  low-quality  boiling 


flux,  high-quality  boiling.  Data  from  Reference  23 


where  utw  is  the  tangential  velocity  at  the  tube  wall.  ufcw  can  be 
calculated  from: 


v  ■  If  r<E™> 


(V-8) 


where : 
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u  la  the  axial  velocity  and  Y  ia  the  number  of  diameters  per 
& 

180-degree  twist  (twist  ratio).  Therefore: 
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Setting  atang  equal  to  g: 
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The  Y  value  calculated  for  the  experiment  with  a  mass-flow  rate 
of  5  X10-3  kg/s  and  an  evaporator  diameter  of  0.25  inches  was  0.66. 
All  attempts  to  make  a  twisted  ribbon  with  this  twist  ratio*  with 
the  facilities  available,  failed.  Therefore,  thiB  method  could  not 
be  experimentally  verified. 


PRESSURE-DROP  CALCULATION 

r  1 

In  a  report  published  by  Bae  et  al.  the  pressure  gradient 
in  a  tube  for  two-phase  flow  is: 


dP 

■ar 


(V-l 2) 


where  dP/dz),  .  „  it>  the  frictional-pressure  gradient,  dP/dz)  is 
rrlct  mom 

the  pressure  gradient  due  to  a  change  in  momentum  as  the  liquid 
changes  to  a  vapor,  and  dP/dz)^rav  is  the  gravitational-pressure 
gradient.  The  frictional-pressure  gradient  can  be  calculated  with: 


4,„ 

■330  Mir) 


.0.0523 

0.105 


X  (l-x)0'117,1 


X  (l-xi0-91*0-86 


(V— 1 3 ) 


The  momentum-pressure  gradient  can  be  found  with: 


VJ  Vft/ 


(V— 14) 


And  the  gravitational-pressure  gradient  can  be  found  with: 


«(oi  - 


(V— 15) 


where  dy  Is  the  differential  vertical  drop. 

These  equations  were  programmed  into  a  computer  routine  and 
compared  with  data  from  the  experimental  apparatus.  The  equations 
predicted  a  good,  average,  pressure  drop.  However,  the  accuracy  of 
the  equations  could  not  be  measured  because  of  the  severe  pressure 
fluctuations  encountered  in  the  test  section. 
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VI  EVAPORATOR  DESIGN  PROGRAM 


The  methods  outlined  in  section  V  were  used  to  develop  a 
program  to  design  straight-tube  evaporators  for  electronics  cooling 
at  10aC  using  either  ammonia  or  freon-1 1  or  at  100°C  using  water. 
Other  fLuids  can  be  used,  but  require  the  addition  of  their  fluid 
properties  to  the  program.  The  user  of  the  program  chooses  a 
working  fluid,  the  total  input  power,  diameter  of  the  tube,  inlet 
temperature,  allowable  maximum  wall  temperature,  and  the  exit 
quality  or  the  mass-flow  rate.  The  program  performs  the  following 
tasks.  It  calculates  the  mass-flow  rote  needed  or  the  exit 
quality;  whether  it  is  possible  to  design  an  evaporator  with  the 
given  inputs;  the  lengths  of  the  one-phase,  subcooled,  and 
two-phase  sections  and  the  total  length;  the  wall -temperature 
differences  for  each  section;  the  average  heat-transfer  coefficient 
for  each  section;  the  pressure  drop  for  each  section  and  the  total 
pressure  drop. 

To  calculate  the  length  of  an  evaporator,  the  following  methods 

can  be  used.  For  an  evaporator  subjected  to  a  constant  heat  flux, 

the  incipient-bolting  point  must  be  calculated  first  using  equation 

V-4.  Then  h  and  h  are  calculated  using  equations  iV-27,  IV-23, 
op  tp 

IV-24  and  V-l,  And  finally,  using  the  allowable  ^T,  average  h^, 

Q,  and  d,  the  length  can  be  found.  For  an  evaporator  subjected  to 
a  constant  wall  temperature,  calculate  the  average  two-phase 
heat-transfer  coefficient  using  equations  IV-23,  IV-24,  and  V-l. 
Since  one-phase  heat-t ransfer  cannot  occur,  this  will  bo  the 
average  heat-transfer  coefficient  for  the  entire,  tube.  Next, 
calculate  the  heat  flux  using  the  given  wall  temperature.  And 
finally,  calculate  the  length.  The  pressure  drop  can  be  calculated 
using  equation  V-l 7. 


Example  1:  Water 

ENTER  FLUID:  WATER < 1 > , AMMON IA ( 2 ) , F-l 1 ( 3  )  1 
ENTER  POWER  <W>  >>1000 

ENTER  DIAMETER  <M>  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K)  >>363.0 
ENTER  DELTA  T  <K>  »10 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  <1>  OR  FLOW  RATE  < 

EXIT  QUALITY  FRACTION  >>.8 


THE  MASS  FLOW  RATE  IS  0 . 540BE-03 < kg/s > 
INITIAL  FLOW  IS  LAMINAR 
BOILING  STARTS  IMMEDIATELY 

ONE-PHASE 

LENGTH  <nt> 

DELTA  T  <  K ) 

h  < W / m ~ 2  K >  233.6  7703. 


SUBCOOL  BOILING 
0 . OOOOE+OO  0.4014E-02 
0 . OOOOE+OO  18.74 


TWO-PHASE 

0.1697 

10.00 

0. 1443E+05 


TOTAL 

0.1737 


DELTA  P  <PSI) 


0. OOOOE+OO  0. 1441E-06 


0 . 3694E-02  0.3694E-0 


Ammonia 


ENTER  FLUID:  WATER < 1 ) t AMMON IA< 2) ,F-1 1 < 3 ) 2 
ENTER  POWER  <W>  >>1000 

ENTER  DIAMETER  <M>  >>1.27E-2 

ENTER  INLET  TEMPERATURE  <k>  >>273 
ENTER  DELTA  T  <K>  >>10 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  < 1 )  OR  FLOW  RATE  <2>?>>1 

EXIT  QUALITY  FRACTION  >>.8 


THE  MASS  FLOW  RATE  IS  0 . 973 1 E -03 ( kg/s > 


INITIAL 

FLOW  IS 

LAMINAR 

ROILING 

STARTS 

IMMEDIATELY 

ONE-PHASE 

SUBCOOL  BOILING 

TWO-PHASE 

TOTAL 

LENGTH 

<m  > 

O.OOOOE+OO 

0. 1626E-01 

0.3433 

0.3596 

DELTA  T 

<  K ) 

O.OOOOE+OO 

10.93 

10.00 

h  <  U/m 

"2  K  > 

100.3 

6377. 

6970. 

DELTA  P 

<PSI) 

O.OOOOE+OO 

0 . 1352E-05 

0 . 1591E-02 

0 . 1593E-02 
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Example  3:  Freon-11 


ENTER  FLUID:  WATER <1 > , AMMON IA< 2 > , F~1 1 < 3 > 3 
ENTER  POWER  <W>  »1000 

ENTER  DIAMETER  <M>  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K)  »273 
ENTER  DELTAT  <K)  “> >  15 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  <1>  OR  FLOW  RATE  (2>?»1 

EXIT  QUALITY  FRACTION  ».8 

THE  MASS  FLOW  RATE  IS  0 . 6377E-02 < kg/s > 

INITIAL  FLOW  IS  LAMINAR 


BOILING  STARTS  IMMEDIATELY 


ONE-PHASE 

SUBCOOL  BOILING 

TWO-PHASE 

TOTAL 

LENGTH  <m> 

0 . OOOOE+OO 

0 . 1967E-0 1 

0.3397 

0.3593 

DELTA  T  <  K ) 

0 .  OOOOE+OO 

17.98 

15.00 

h  (  W  /  to  ”  2  K) 

35.73 

3880. 

4650. 

DELTA  P  <PSI) 

0. OOOOE  +  OO 

0.9424E-05 

0.8642E-01 

0 . B643E-01 

Example  4 :  Subcooled  We ter 


ENTER  FLUID:  WATER < 1 ) , AMMON IA< 2 ) , F-l 1 < 3 > 1 
ENTER  POWER  <W)  >>500 

ENTER  DIAMETER  <M)  »1.27E-2 

ENTER  INLET  TEMPERATURE  <K>  »330 
ENTER  DELTA  T  <K)  >>15 

DO  YOU  WANT  TO  ENTER  EXIT  QUALITY  (1)  OR  FLOW  RATE  <2>? 

ENTER  MASS  FLOW  RATE  <k.3/s)  >>6.0E-3 

ENTIRE  SECTION  IS  SUBOOOLEP 


VII  CONCLUSIONS 


1.  From  the  review  of  the  forced-convective  boiling  process: 

a.  The  gravitational  acceleration  has  very  little  effect  in 
the  subcooled-boiling  and  two-phase,  f orced-covective  regions. 
Therefore,  the  lack  of  a  gravitational  acceleration  (o-g)  should 
not  affect  the  calculation  of  local  heat-transfer  coefficients. 

b.  The  gravitational  acceleration  does  affect  the  saturated 
nucleate-boiling  region.  However,  since  the  gravitational  effect 
tends  to  balance  out  around  the  circumference  of  the  tube  and  since 
this  region  is  generally  very  short,  the  use  of  equations  developed 
for  l-g  accelerations  can  be  used  to  calculate  the  local, 
heat-transfer  coefficient  without  much  loss  of  accuracy. 

2.  From  the  flow-visualization  study,  placing  a  twisted-tape 
insert  into  the  evaporator  should  stabilize  tin  flow  by  inducing  a 
large  centrifugal  acceleration  which  prevents  slug  flow  from 
eccuring.  The  twist  ratio  must  be  large  enough  to  ensure 
stratified  flow  under  all  conditions.  This  twist  ratio  can  be 
calculated  using  equation  V-ll. 

3.  The  length  of  a  straight-tube  evaporator  can  be  calculated  for 
a  constant  heat  flux  or  a  constant  wall  temperature  by  using  the 
methods  developed  Ln  this  study  and  summarized  in  Section  VI. 
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APPENDIX  A 


This  appendix  shows  how  the  calculations  were  done  to  obtain 
the  solid  lines  on  the  graphs.  More  specifically,  sample 
calculations  are  shown  for  Figures  8,  9,  and  10.  In  addition,  a 
sample  twist  ratio  calculation  is  done. 


«r  j  «r „  a  -  t  *, « 


/-'/-V'  *•>  *■'  /  ■*  \  \\\ ,  -*>  r 


q  «  1.06X106  W/m2 


SAMPLE  CALCULATION;  Subcooled  Boilings  Figure  9 


One-Phase  Forced  Convection: 
Data:  VTFL  *60  K 

Tfl  *  30  C 
u  -  0.5  m/s 
d  *  1.2X10"2  m 


«  963.3  kg/m3 
p  =  1*. 71X10“^  kg/m  s 
k  *  O.65I+  W/m  K 
Pr  ■  3.01 


Equation: 


h 


FC 


0.23Re°'8Pr0,1< 


Calculation: 


ref 


w 


T  +t_t 
-  FL 


90+30 

2 


60  C 


nr  pud  _(983.3kg/m3)(0.5m/s)(l.2X10'2m) 
p  1*.  71X10  Ag/m  s 


Re  =  1.253X10** 


hFC  *  0.023(1. 253X10U)0,8(3.01)°-14 


0.651+W/mK 

1.2XKT2m 


hpc  -  3697  W/m2  K 

q  -  LfcAT 


q  a  (3697W/m2K)(60K)  *  2.22X105  W/m2 


Incipient  Boiling  Point: 


Data: 


■  0.68  w/ffi  K  0  -  5.89X10“£  x/ro 

p.  s  953.0  kg/m3  Pr  s  1.66 
hf  =  2.258X10  J/kg  Vfg  -  1.6719  sT/kg 
U  -  2.8X10"\g/m  s  A  *  1.2X10“2  m 


Equati on : 


,  Vrc 


3St’IB  *,t>, 

H  fg 


Calculation: 


He  =  (958kg/ni  ) (0 .5m/s ) ( 1.2X10"  m) 

y  2. 8X10“  Kg/a  s 


Re  »  2.053X10 


0  8  Oik 
hpc  *  0.023ReU-°Pr-"  §• 


h__  *  0. 023(2. 053X10 11  )0,8(  1.66 )°- °-68w/n^ 


1.2X10  n 


h  «  1198  W/ra  K 

r  C 


sat , IB 


=  (8)(5.89X10'~2N/m)(373K)(l.6719in3/kR) 
(0.68V/BK) (2. 258X10^ J/kg) 


X  (Ul7W/m2K)(l.66)2 


ATsat,IB=2'37K 


Subcooled  Boiling: 


Data: 


py  -  0.6  kg/ra3  P  =  2X105N/m2 

Cp  =  1.21X10 3  J/kg  K  AT  =  20  K 


wn 


V.VAA 
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Equation: 


$ 

I 

/■ 

$ 


o 


h  *  0 
NB 


ft-) 


0.31 


4fg* 
0.31 


hNB  S  °-25 


•  _g  .0,69 

a  fiBtf/mK  a (1.2X10  »)_ - j-— - 


/  0.68W/mK  a(  1.2X10  - - - 

\1.2X10“2m  (2.258X10  J/kg) (2. 8X10  kg/Bs)/ 

x  /(4.21X103J/kftK)(2.8Xl0^kg/ms)\0,69 
\  o.68VJ/»X  / 

( ( 2XlO^/m2 Hi. 2X10"2m )N° * 31  /  ?58k£.^  ,  A 
X  V  5.89X10“2N/m  /  \0.6kg/m  / 


0.69 


hNB  =  3.02q 


9  "  hNBATsat 


q  -  3.02q°'696Tsat 


9*  35.4ATsa3-23 


q  *  35.14(20  K)3,23  -  5.65X105  W/m2 


s  ' 

sample  CALCULTION:  Forced-Convective 

Boiling 

« 

Data:  k  -  0.68  W/m  K 

hfg " 

~-r 

p^  -  958  kg/a3 

Op  * 

i 

p  »  0.6  kg/m3 

0  s  ! 

jy 

'i 

-  2.8X10"1*  kg/m 

s  U  * 

LO3 

v-2. 
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B-fc.  B  UM  u, V- 1 .'S  TV  .  L-Ts-t r* 4 r* ,  »*-.-aA.-*w»  .  -«  », ,  j 


d  si .91X10  m 


q  *  1*3525  •:/-/ 


Equations: 


Calculations: 


P  ®  1.1*5X105  N/rn2  Pr  *  1.66 


vf  “1.6?  m3/kg 


x  *  0.03 


,  .  v0.69/ojS0.31  , 

k/qd  \  /Cp)j\  /  Pd\  /p  \ 

=  °-25Kv)  C~)  (.-)  fo-v 


,  8pT^i-V?.q  (Pr)2 
sat,IB  Vfg 


hpc  =  0.023Re°,8Pr°-l<^' 


_  ■.,tVvti> t  • 

qFC  Fgd  ATsat 


(one-phase ) 


(two- phase) 


V  •..(*- fe“)) 

.-v  -6^)') 


Re  =  md  =  ( ,21kg /s  )  (1.91X10"  tn) _ 

Ay  7T/l*(l.9lX10"2m)2(2.8X10“\g/ms) 


Re  s  5X10 


h  =  0 . 02  3 ( 3 X10 ^  )°  ’ 8  ( 1 . 66 )°  ‘  ^ 

1.91X10“  m 


=  5760  W/tn  K 


(aW5.89X10"2M/m)(373K)(i.67nr/kgj 


sat, IE  (0.68W/mK)(2.758X10®J/kg) 


X  (5760W/a2K)(l.66)2  •  2.71  K 


MB 


.  0.25/kli£^WM) _ nr - ) 

\(  2 . 25 8X10® J/kg ) ( 2 . 8X10"4kg/ms )/ 


0.69 


( It .  ?lXlQ3J/kgK)  ( 2 , 8X10~I*kg/as  )^  0,69 


0 ,60W/mK 


h  liSY103N/m2)(l.91X10"^m) 


\  5. 


-2  AO.  31 


.  89X10“  N/m 

X  /iMSs/al.  1V'31 

\0.6kg/nr  / 


9 


NB 


2.735q°‘69  *  2.735d3525W/a2)0,69 


iHB  «  **3Ul  U/mcK 


/0.6kg/m3>\0,Y  2.8XloAg/ms  \0,1/d-Od 

'\958kg/o3/  \1. 27X10“ ^tg/rns/  \  0.62 


Xt  \958kg/o3 


X  =  0.7788 

X»  u 


Re  =  R.(l  -  *)  *  *>.85X105 


F(XM)  •  0.15 


c- 


7788 


+  2.0 


jlJ{ 


0.32 


.7788; 


F(Xtt)  *  0.5176 


F  =  5(1.66)  +  51n(l  +  5(1.66))  +  2.51n(0.0031 


X  (M5X101*)0*812) 


F2  =  26.91 


*FC 


„  (**.  85X10^  )0, 9(  1, 66 )(  .51T6)(.068w/mK) 

( 26 . 9 1 ) ( 1 . 9 lX10"2m ) 


AT 


sat 


qrc  -  lfW3ALat 


q  =  U3525  =  187I1MT  -f  U3UlATsat 


e 


2.71 


AT 


sat' 


AT 


sat 


2.1*9  K 


since  ATsat  ^  *  2.71  >  ^TEat»  boiling  is  suppresed  and  q  *  qFC> 
There  fore: 

h  *  1371*3  W/m2K 


■t 


SAMPLE  CALCULATION:  Twist  Ratio 

Data:  A  -  5X10“  3  kg/s  d  =  0.25  in  =  6.35X10“3m 

=  955  kg/m3  0  =  5.89X10“2  N/m 


Equation: 


Calculation: 

u  =  Ji  s  5X10~3kg/s _ 

a  pA  (955kg/m3)(1T/l*(6.35XlO"3m)2) 

u  =  1.65X10-1  m/s 
a 
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Y  “  0,66  diameters/180  degree  twist 
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